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Abstract 

We investigate the effects of the gravitational field on the quantum 
dynamics of non-relativistic particles. We consider N non-relativistic 
particles, interacting with the linearized gravitational field. Using the 
Feynman - Vernon influence functional technique, we trace out the 
graviton field, to obtain a master equation for the system of particles 
to first order in G. The effective interaction between the particles, 
as well as the self-interaction is non-local in time and in general non- 
markovian. We show that the gravitational self-interaction cannot 
be held responsible for decoherence of microscopic particles due to 
the fast vanishing of the diffusion function. For macroscopic particles 
though, it leads to diagonalization to the energy eigenstate basis, a 
desirable feature in gravity induced collapse models. We finally com- 
ment on possible applications. 
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1 Introduction 



There has been recently a considerable interest in the application of the 
influence functional technique in the study of non-equilibrium systems in 
physics. Besides quantum Brownian motion 0, 0, ^, ^] for which the method 
was initially developed, it has been applied to the modelling of particle- 
field interactions |0, radiation damping [0], black-body radiation and 
most recently to non-inertial detectors coupled to a scalar field 0. It is 
one of the most powerful techniques to obtain master equations, when the 
coarse-graining comes from the splitting of degrees of freedom to system and 
environment. 

In this paper we apply the technique in another system of N non- 

relativistic particles coupled to linearized gravity. A motivation for this is 
the possibility that gravity induces decoherence on the particles' states. This 
is a suggestion made on different contexts on fundamental irreversibility in 
quantum mechanics |]Tl], 0, |T2|. The weakness of the coupling suggests that 
the probable decoherence time should be very large, but the particular form 
of the coupling (quadratic to momentum) and the possibility of persistent 
noise might give rise to observable consequences. 

In addition, the model we present here can be generalized in a straightfor- 
ward way to obtain a description of systems of quantum mechanical detectors 
of gravitational waves. 

Our model consists of non-relativistic particles coupled to the linearized 
gravitational field, which is assumed to be initially at its vacuum state. We 
argue that a factorizing initial condition is, in contrast to quantum Brown- 
ian motion, well suited for our system. The modes of the graviton field are 
bounded in energy by an ultraviolet cut-off A, which on physical grounds 
should be much smaller than the Compton wavelength of the particles. In 
addition, we assume that the particles are almost stationary. Our analy- 
sis resembles, in a way, the one of 0. Like them, we obtain correlation 
kernels describing a non-local interaction between the particles. The influ- 
ence functional we construct is rather different from the ones considered in 
the literature, due to the particular features of the gravitational coupling 
which is quadratically coupled to momenta. The result of our analysis is the 
non-markovian master equation (3.9) . For the case of a single particle, it is 
simplified significantly (4.1). We see that the dissipation and diffusion are de- 
termined solely from the Hamiltonian operator. We can interpret our results 
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as a continuous monitoring of the energy of the particle by the gravitational 
environment. 

The diffusion function, which is responsible for decoherence vanishes at 
long times and it turns out, that unless we consider macroscopically massive 
bodies, the rate of gravitationally induced decoherence is extremely small. 
This is a desirable result in connection with the gravity-induced collapse 
models. 



2 The model 

Consider N particles on a 3+1 dimensional spacetime, moving on trajectories 
(x„(r„), t„(r„)) parametrized by the proper times r„ so that t„(r„) is a strictly 
increasing function of r„ 0. We assume that the gravitational interaction is 
very weak and therefore work in the linearized approximation. That is, the 
metric is: 

9t^u = V/iiy + hf^u (2.1) 

with rj^i, the Minkowski space metric. 

We take the non-relativistic limit for the particles, that is, we assume 
that there exists a frame, with respect to which they are almost stationary 
and therefore can write their trajectories as (a*^^-) + x^(^^^{t),t), having identi- 
fied the global time coordinate t with the proper-time of the particles. We 
assume that | X(„) | is much smaller than the distance between any two par- 
ticles dnm = |cL(n) — 8L(^m)\ ■ This is a good approximation as long as dnm is 
much larger than the maximum wavelength of the graviton field that can be 
excited. Essentially, we consider the particles moving around some fixed sites 
coordinatized by a(„), so that their individual motion does not significantly 
change their distances. In any case, this approximation does not affect at all 
the discussion on the self interaction of the particles through the gravitational 
field. 

We work in in the transverse-traceless gauge for the linearized gravita- 
tional field (/lo/i = 0, h^i,= 0, h] = 0). Under these approximations , the 
total action of the system for evolution from global time t = to t = T, 
reads 

Slot = Sgr + Spar + ^int (2-2) 
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where 

Sgr = ^fdtj d^xh^,,ph^^'' ^^J^dtJ d^xikjh'^ - h^j,kh'''') 

(2.3) 

Spar = J2f dt:^SiAn)iU (2-4) 
„ JO Z 

n "'O 

Note that we have set h = c = m = 1. 

We expand the graviton field in normal modes: 

^'^^""^ = / (0 ^(^Sc^^kx + g(^.^r^kx)4?, (2.6) 

The polarization matrices A^^j ( r = 1,2) are traceless and transverse 
and can be chosen to satisfy : 

4-X'^' = 5..'(^.^-^) (2.7) 

E4-i4;\ = (% - ^)(5.). - ^) := T-'='(k) (2.8) 
The gravity part of the action therefore reads: 

= 2^ /o7 (0 + '^'^-'^ + + ^2.9) 

This is just the action for two massless scalar fields propagating on Minkowski 
spacetime. 

We now write the coupling part of the action 

^llo'^'Iwr^^ E(?Sc««ka(„) + q^^sink^^^))4lx^x^ (2.10) 

where within our approximations we ignored the X(„) terms in the trigono- 
metric functions. 
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By using the collective index a to include the k,r and the indexing of our 
oscillator by 1 or 2, we write : 

Sgr + S.nt = £T.[^iQl + ^ya)+^M (2-11) 

where 

jS = coska^n)A^^lx'x^ (2.12) 

4^ = sznka(„)Agi;^i;^' (2.13) 

and cjk = |k|. This is just the action of a collection of forced harmonic oscilla- 
tors. Therefore the total action is that of a collection of N non-relativistic free 
particles interacting with a bath of harmonic oscillators, through couplings 
depending quadratically on the velocity. 

The tracing out of the graviton modes can be done exactly since the 
path integral is a gaussian with respect to them. We compute the influence 
functional: 

^[x(t),x'(0]= (2.14) 
Ujdqfj dqf I dq^ j dq^5{q^ - | Dq^{t)Dq'p{t') 

e^-p[iSgr[qa{t)] + iSint[qa{t),:!^{t)] - iSgr[(^^{t')] - iSint[q'a{i')-,^{t')]] 
Po(go,xo,g^f, /6/x'o) 

where the integration is over the paths satisfying: g"(0) = , q'^iT) = qj, 
g'"(0) = ^q", q°'{T) = q'j'. Here po is the density matrix of the total system. 
The path integrations can be carried out exactly, to obtain: 

^[x(t),x'(t')] = 

H{T) exp[- f ds f ds\J^ + j;)(s) smu^{s - s'){J^ - J'J{s') 

^ ZCUa Jo Jo 

f ds fds'iJ^ - J'Jis) cosa;,(s - s'){J^ - J'J{s')] (2.15) 

^ ZUJa Jo Jo 

In deriving this we have assumed that at t = the states of the particles 
and of the graviton field were uncorrelated and that the field were on its 
vacuum state, i.e. 

*[M=CexpE-^q^] (2.16) 
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This initial condition is usually considered unphysical in quantum Brownian 
motion models. We believe that it is actually a quite good one for the case 
of gravity. Graviton modes are excited only by non-stationary particles. 
Therefore, this initial condition reflects an operation on the particles of a 
very fast acceleration just before t = 0. 

Substituting the expressions for the currents Ja into the influence func- 
tional we get: 

^[x,x'] = iV(T)exp[z^ f^ds f ds\x\^)xl^ + x\^)X'l^){s) (2.17) 

7(n)(m)('^ ~ )i^tm)^[m) ~ ^'(m)^'(m) ) (■^ ) 
-J2 [ / ^^'(i(n)i(n) -^'{n)a;'(n))(«) 

The kernels 7(n){m) and rj(n){m) are given by the expressions: 

^(5U(^) = ^/ ^sin|k|.cosk(a„-a^)r^'='(k) (2.18) 
C r d^h 

= ^2 J ^cos|k|.cosk(a„-a^)r^'='(k) (2.19) 

These are the dissipation and noise kernels, similar to the ones derived in 
[y] for the case of detectors minimally coupled to a scalar field. For n ^ m 
they describe the dissipation and diffusion induced on the particle n from the 
particle m , while for n = m they contain the effects of the self-interaction 
of the particle through its interaction with the gravitational field. 

In order to keep them finite, we have to restrict the integration range to 
values of |k| smaller than a cut-off A. This is natural, since we do not expect 
the non-relativistic particles to excite graviton modes with arbitrarily high 
energy. In fact A should be much smaller than the Compton wavelength of 
the particle. This is in accordance with our previous approximations, since 
the distance between any particles remains much larger than their Compton 
wavelength. 
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In the particular case n = m we can perform the angular integrations in 
spherical coordinates in the equations for the kernels and obtain: 

^Sin)i') = ^^'''' C^kksmhs (2.20) 

€)Us) = ^5^''' j^dkkcosks (2.21) 

We note that by taking the cut-off to infinity, the dissipation kernel be- 
comes essentially the derivative of a delta-function, as in the well studied 
case of quantum Brownian motion with ohmic environment. The correspond- 
ing semi-classical equations for t >> A^^ can be found using the standard 
procedure 0, ||, 0: 

2C 

x' + —5,,kix'x^x' = + iH'')iki (2.22) 

with ^kiif) a stochastic force determined by the correlator: 

{Ut)iki(t')) = r,'''\t - t') (2.23) 



3 The master equation 

Having obtained an expression for the influence functional we can compute 
the reduced density matrix propagator: 

J(xf,Xf,t|xo,Xo,0) = y" j L)xi:)x'exp(iS'par[x] - i5'par[x'])J'[x,x'] (3.1) 

where the integration is over all paths x(s), x'(s') satisfying: x(0) = xq, 

X'(0) = X;,, X(t) = Xf, X'(t) = Xf. 

The knowledge of the reduced density matrix propagator enables us to 
construct a master equation. Our system is characterized from the non-local 
dissipation and diffusion in the influence functional, and the coupling which 
is quadratic to the velocities. Because of the peculiarities of the latter, the 
general method of Hu,Paz and Zhang is not applicable here. Instead we 
compute the influence functional perturbatively (first order in G) and use 
the Feynman prescription for the determination of the master equation. 
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Our starting point is the density matrix propagator for the free particle 
under external forces F(s), F'(s): 

J(°)[F, F'](xf, x'f, t|xo, x'o, 0) = ^ exp[^(xf - Xo)^ - ^(x'f - x^,)^ (3.2) 

+ -xo / dssF(s) - -x'o / dssF'(s) 
t Jo t Jo 

+-Xf / ds{t - s)F(s) - -Xf / ds{t - s)F'(s) 
t Jo t Jo 

+ - f ds r ds's'{t-s)F{s)F{s')-- f ds r ds's'{t-s)F'{s)F\s')] 
t Jo Jo t Jo Jo 

The perturbation expansion of the propagator is writen then formally: 

J(xf,Xf,t|xo,Xo,0) = (3.3) 



To first order in G we obtain: 

tlXn.xLO) = ^ 



c 

J (xf , x^, t |xo, x;,, 0) = - exp J2 [^C^ii^iC/Jn^U (2-^) 



(n)(m) 

~ ^o)(n)(2^/ - a;o)(„)5mn - -^^ijix'f - a;o)(„)(a;j - a;o)(^)5„ 

-y (3/ - 4ig)g}^)5i^(x/ - Xo)5„)(x/ - Xo)J„)5„ 

G 



— {x'f — a;o)(^)(a;y — a;Q)-^^)(a;y — a;o)('^)(a;j — a;Q)(^)] 



G 

{xf — Xo)(„)(a;j — Xq)-'(^^-^{x ^ — Xq)(^^)(xj — Xq)!^) 
-2(a;/ - Xo)(„)(a;j - Xo)[^^{x'f - x'Q)\^){x'f - x[,)[„)]] 
where / and g are functions of time: 
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and in particular: 

-i^^'^'r*^ (3.8) 

The standard prescription for the derivation of the master equation from 
the reduced density master propagator consists of taking its time derivative 
and using identities relating xq and Xq with the action of derivatives with 
respect to Xf and Xj . For the interested reader, we hst the relevant identities 
in the appendix. 

After some calculations, the master equation turns out to be (inserting 
back h,m and c): 

|,.X:|.(i_5™,„,(*))(^-a|^)/> (3.9) 

r)'^ 

R^jkl (t)( I 

-2 ^ )p 

This is the main result of this paper: the master equation for N non- 
relativistic particles interacting through linearized gravity. The gravitational 
field induces a renormalization in the mass of the particles, modifies the 
dynamics so that they become dissipative and is responsible for noise. These 
three effects are contained in the functions 5m{t), a{t) and (5{t) respectively: 

= ^^J^5(n)^ij-^w (3.10) 



a 



ijkl (,\ fijkl J.2 



(n)(m)(^) - ^^5f{n){m)^ (^-H) 
AC 1 

/^;5U=^(^^+2^^')hU (3-12) 
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4 One particle 



An interesting case is that of a single particle. Since then the functions in the 
master equation are totally symmetric in the spatial indices, we can, without 
loss of generality, consider it constrained to move in only one dimension. The 
master equation reads then in operator form 

-p = -'-[Hji, p] - ia{t)[H% p] - (3{t)[Hn, [H^, p]] (4.1) 

and depends explicitly only on the renormalized Hamiltonian H^. We can 
verify that this form of master equation (in particular the noise part) is 
particular to the free particle case. For an harmonic oscillator we would get 
an extra dissipation and diffusion term due to the coupling of the particle's 
position to the graviton oscillator Hamiltonian, and of form similar to the 
one derived in for quadratic coupling to position. 

The diffusion coefficient f5{t) exhibits a "jolt" for times of the order of 
A~^. In quantum Brownian models this is a cause of rapid decoherence of 
the density matrix of the particle, and diagonalization in a basis determined 
by the coupling to the environment. Our particular form of the diffusion 
terms tempts us to propose that it should lead to diagonalization of the 
particle's density matrix in the energy eigenstate basis. But we have to take 
into account, that the coupling is extremely weak and that after the jolt the 
diffusion coefficient falls to zero, quite slowly actually since it goes at most 
like 1/t. 

We can give an estimation of the decoherence in the energy by approx- 
imating l3{t) with a constant of the order of for times of the order of 
A^^ and zero afterwards. We borrow some ideas from the quantum state 
diffusion picture of quantum mechanics |T^, n, |T3|. At the times that f3{t) 



is constant, we have a unique unravelling of the density matrix into states 
evolving stochastically in a Hilbert space. It is straightforward to show |]13| 



that an initial wavepacket with energy spread AEq will emerge after the jolt 
with spread AE given by: 

' (4.2) 



{AEf (AEo)2 nc> 

lass m a good upper bo 
classical gravitational self energy of a mass distribution localized within the 



For a single particle of mass m a good upper bound on UK is ^^^-^ : the 
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Compton wavelength of the particle. This means that: 



(AE)2 (AEo)2 ~¥^- 



-2-1^ (4.3) 



This is an extremely small quantity, when considering microscopic particles 
(even in atomic scales). On the other hand, for macroscopic and even meso- 
scopic particles the right hand side is quite large and we expect a loclalization 
of the particle in its energy eigenstates. For instance a particle with mass 
m = 10~^gr and irrespectively of its initial configuration, will emerge after 
localized in an energy eigenstate with spread of the order of O.lMeV, 
which is a tiny portion of its kinetic energy. But in this case, the gravity 
induced decoherence is in general, hidden beneath the effects of other types 
of environment |jT6[. In any case, this result is in good agreement with the 
assumptions of the gravitationally inducced collapse models. 

These features were , more or less expected, since gravity couples very 
weakly and its strength increases with the mass of the interacting bodies. 
Still,there was the possibility, that a persistent noise source might induce de- 
coherence even in microscopic systems, despite the weakness of the coupling. 
Note, that our analysis based on the linearized approximation, does not rule 



out the possibility that highly non-linear Planck scale processes [|lT| might 
be a source of noise, giving rise to decoherence at smalles mass scales. 

The dissipation function a{t) approaches asymptotically a constant of the 
order of The overall picture we get, is that of a particle continuously 
dissipating energy and suffering at early times noise from the environment 
until it becomes correlated with the gravitational field. 



5 Conclusions 

We have studied the quantum theory of N non-relativistic particles, coupled 
to the linearized gravitational field using the influence functional formalism. 
Our main result was the master equation (3.9) containing information of 
non-local interaction between the particles. We should note that the grav- 
itational field, being coupled quadratically to the velocities gives a rather 
unusual expression for the influence functional. This results in a master 
equation, where both dissipation and diffusion are determined uniquely by 
the Hamiltonian operator. This is in accordance with our intuitive feeling. 
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that the gravitational field acts like continuously "measuring" a particle's 
energy. 

One of our motivations for this work, was to establish whether we can 
consider the gravitational field as a source of fundamental decoherence in 
quantum mechanics. The answer comes out negative for microscopic systems, 
but systems with large mass seem to decohere within a fast rate, in the energy 
eigenstate basis. In addition, it might be interesting to examine, the evolution 
of a single particle under the action of a particular matter distribution. The 
formalism we used can be extended with slight modifications to cover this 
case. We can, for instance, consider almost stationary cosmic dust and even 
a cosmological spacctime. The collective effect of matter plus gravity might 
give a strongest decoherence to the particle. 

In addition, it would be of interest to study the response of a system of 
detectors to different initial conditions for the graviton field. The case where 
a number of modes is excited, seems very interesting. The information of the 
state of the field should be encoded in the correlation kernels of the particles, 
from the time evolution of which we would be able to determine the presence 
of the graviton fields. This might give a nice toy model for detectors of 
gravitational waves. 
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7 Appendix 

We give here the identities that enable us to compute perturbatively the 
master equation. We give the form for the case of one dimension and one 
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particle. The generalization is straightforward. 

{xf - xo)'^J{xf, x'f, t I xo, x'o, 0) = (7.1) 
[-t\l - 2zG(9/ - - ^) - 2Gt^(^ + ^) - 2Gt(3/ - 4tg) 

+ ^(/ - - a^o)^ + -J-9{.Xf - xofi^x'^ - x'Qf]J{xf,Xf,t \ Xq,x'^^,Q) 

{xf - xo)^ J(x/, x'f, t I Xo, Xq, 0) = (7.2) 
~ ^'^'94 ~ 3iV(^/, 4, 1 1 Xo, 4, + 0{G) 



{xf - Xo) {x'f - x'o) J {xf, x'f, t I xq, x'o, 0) = (7.3) 

^^'dx]d^^ + '^'dx^f ~ '^'d^ + ^ ' °^ + ^^^^ 

We should keep in mind that eventually, we keep terms to first order in G. 
The expressions for the primed quantities are obtained by permutation of 
primed with unprimed ones and complex conjugation. 
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